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We report the observation of decays B° -^ D),*''^ti~ and B° -^ Ds~K'^ in a sample of 230 x 10'^ 
T(45') — > BB events recorded with the BABAR detector at the PEP-II asymmetric-energy e'^e~ 
storage ring. We measure the branching fractions 3{B'^ — > Dfir') = (1.3 ±0.3 (stat) ±0.2 (syst)) x 
10-^ B{B° ^ D7K+) = (2.5±0.4 (stat) ±0.4 (syst)) x 10-^ B{B° -^ Dt+n-) = (2.8±0.6 (stat) ± 
0.5 (syst)) x 10-^ and BiB" -^ D*' K+) = (2.0 ± 0.5 (stat) ±0.4 (syst)) x 10~^ The significance 
of the measurements to differ from zero are 5, 9, 6, and 5 standard deviations, respectively. 



PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



Within the Cabibbo-Kobayashi-Maskawa (CKM) 
model of quark- flavor mixing ,1] , CP violation manifests 
itself as a non-zero area of the unitarity triangle |2| • One 
of the important experimental tests of the model is the 
determination of the angle 7 = a.rg{-Vu,dV*f^/VcdV*,^) of 
the unitarity triangle. A measurement of sin(2/? ± 7) 
can be obtained from the study of the time dependence 
of the B^,B'^^D^*^~tt'^ J2I decay rates, and specifically 
of the interference between the CKM-favored B° de- 
cay amplitude and CKM-suppressed B^ amplitude [J|. 
The first measurements of the CP asymmetry in decays 
i3"^D'^*^±7r* have recently been published 5]. 

The measurement of sin(2/3±7) in i3°^£'(*)±7r^ de- 
cays requires knowledge of the ratios of the decay ampli- 
tudes, r(i:i(*)7r) = \A{B'^~^D^*'>+Tr~)/AiB°^D^*'>-TT+)\. 
However, direct measurement of the branching fractions 
B{B*^^D'^*'>^n^) is not possible with the currently avail- 
able data sample due to the presence of the overwhelm- 
ing background from B°^-_D*-*^+7r~. However, assuming 
SU(3) fiavor symmetry, r{D^*^Tr) can be related to the 
branching fraction (BF) of the decay B'^ 



Di*^+7r- 



i: 



r{D'^*K) ^tane*. 



foM B{B"^D. 



W + TT-l 



/„(., V B{B" ^ D(*)-7T+) ' 
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where dc is the Cabibbo angle, and /£>(*) //r,(«-) is the 

ratio of I?^*' and Dg* meson decay constants [Sj. Other 
SU(3)-breakingeffects are believed to affect r(D^*^Tr) by 
less than 30% 0. 

Since B^ —>■ Ds 7r~ has four different quark flavors 
in the final state, only a single amplitude contributes 
to the decay (Fig. QJ;). On the other hand, there 
are two diagrams contributing to B^ -^ D''*'^tt^ and 
S" -^ Z?^*'+7r~: tree amplitudes (Fig. ^,b) and color- 
suppressed direct W^-exchange amplitudes (Fig. ^ije). 
The latter are assumed to be negligibly small in Eq. ^. 
The decays B° -^ D'f^~ K+ (Fig. [Tf) probe the size of 
the H^-exchange amplitudes relative to the dominant pro- 



cesses B° -> D'^*'>-'k+. The rate of B^ 



D 



(*)- 



K+ de- 



cays could be enhanced by final state rescattering 7j , in 
addition to the PF-exchange amplitude. 

The branching fractions B{B^ -^ D^tt^) and 
B{B° — > D~K~^) have been measured previously 
by the BABAR [3 and Belle P collaborations, but the 
decays B" -^ D*^tt~ and -B" -^ D*^K~^ have never been 
observed. In this Letter we present new measurements 




FIG. 1: Dominant Feynman diagrams for (a) CKM-favored 
decay -B° -^ D^*'~n^, (b) doubly CKM-suppressed decay 
B" — > D'*'^7r~, and (c) the SU(3) flavor symmetry related 
decays B^ -^ Ds* tt"; (d) the color-suppressed VF-exchange 



contributions to _B" - 
decay B° -^ oi*^' K"' 



D<*'-7r+, (e) B° -^ D^'^+n-, and (f) 



iW- 



of the decays B° ~* D)*>^tt- and B" -^ DT'~ K+ . The 
analysis uses a sample of 230 x 10^ T(4S') decays into 
BB pairs collected with the BABAR detector at the 
PEP-H asymmetric-energy B factory |lX| . 

Since the BABAR detector is described in detail else- 
where ;11|. only the components that are crucial to this 
analysis are summarized here. Charged particle track- 
ing is provided by a five-layer double-sided silicon ver- 
tex tracker (SVT) and a 40-layer drift chamber (DCH). 
Ionization energy loss {dE/dx) in the DCH and SVT 
and Cherenkov radiation detected in a ring-imaging de- 
vice are used for charged-particle identification. Photons 
are identified and measured using the electromagnetic 
calorimeter (EMC), which is comprised of 6580 thallium- 
doped Csl crystals. These systems are mounted inside a 
1.5 T solenoidal superconducting magnet. We use the 
GEANT4 [121 software to simulate interactions of parti- 
cles traversing the BABAR detector, taking into account 
the varying detector conditions and beam backgrounds. 

We pre-select events which have a minimum of four 
reconstructed charged tracks and a total measured en- 
ergy greater than 4.5 GeV, determined using all charged 
tracks and neutral clusters with energy above 30 MeV. In 
order to reduce "continuum" e^e^ ~* 99 (9 = u, d, s, c) 
background, the ratio of the second and zeroth order Fox- 
Wolfram moments [l3| must be less than 0.5. 

Candidates for D+ mesons are reconstructed in the 
modes D+ -^ (I)tt+J{°K+ and K*°K+, with (j)-^K+K~, 
/-CS-^TT+TT-, and K*°'^K-iT+ 



The K2: candidates are 



reconstructed from two oppositely-charged tracks, and 
their momenta are required to make an angle Inflight | < 
11° with the line connecting their vertex and e+e~ in- 
teraction point. All other tracks are required to origi- 
nate from the e+e~ interaction region. In order to reject 
background from D+^iir°7r+ or ^*"7r+, the K+ can- 
didate in the reconstruction of Df^K'^K+ or K*°K+ 
is required to satisfy positive kaon identification criteria 
with an efficiency of 85% and 5% pion misidentification 
probability. The same selection is used to identify kaon 



daughters of the B mesons in decays S —f D 



(*)- 



K+. In 



all other cases, kaons are not positively identified, but in- 
stead candidates passing pion selection are rejected. Such 
"pion veto" has an efficiency of 95% for kaons and 20% 
for pions. Pion daughters of B mesons in the decays 
B" ^ Ds TT~ are required to be positively identified. 
Decay products of <j), K*^ , K'^, Df, and i?° candidates 
are constrained to originate from a single vertex. 

We reconstruct D*^ candidates in the mode 
D*~^^Df'-f by combining Df and photon candidates. 
Photon candidates are required to be consistent with an 
electromagnetic shower in the EMC, and have an energy 
greater than 100 MeV in the laboratory frame. When 
forming a D*~^ , the D'^ candidate is required to have in- 
variant mass within 10 MeV/c^ of the nominal value [lj|. 

After an initial pre-selection, we identify candidates 
for 5° -> Di*'^^7r- and B° -> D^*^' K+ using a likeli- 
hood ratio Rl — £sig/('Csig -I- -Cbkg), where Csig ~ 
rii 'Paig{xi) is the multivariate likelihood for signal events 
and £bkg — Ili ^bkg(a^i) is the likelihood for background 
events. The ratio R^ has a maximum at i?^ = 1 for 
signal events, and at i?^ = for background originat- 
ing from continuum events. It also discriminates well 
against generic B decays without a real D'^ meson in the 
final state. The likelihoods for signal and background 
events are computed as a product of the probability den- 
sity functions (PDFs) Vsig{xi) and Pbkg(a;i) for a num- 
ber of selection variables Xi'. invariant masses of the (/>, 
K*^ and K^ candidates, y^ confidence level of the ver- 
tex fit for the i?° and D'^ mesons, the helicity angles 
of the (f), K*°, and D*~^ meson decays, the mass differ- 
ence Am(D*+) = m{D*+) - ■m{D+), the polar angle 9b 
of the B candidate momentum vector with respect to 
the beam axis in the e"'"e~ center-of-mass (cm.) frame, 
the angle 9t between the thrust axis of the B candidate 
and the thrust axis of all other particles in the event in 
cm. frame, and event topology variable J-. Correlations 
among these variables are small. The helicity angle 9h is 
defined as the angle between one of the decay products 
of a vector meson and the flight direction of its parent 
particle, in the meson's rest frame. Polarization of the 
vector mesons in the signal decays causes their helicity 
angles to be distributed as cos^ 9h {4> and K*°) or sin^ 9h 
{D*~^), while the random background combinations tend 
to produce a more uniform distribution in cos 9h- 



Variables cos 9b, cos 9t, and JF discriminate between 
spherically-symmetric BB events and jetty continuum 
background using event topology. BB pairs form a nearly 
uniform |cos^y| distribution, while |cos^y| distribution 
for the continuum peaks at 1. A linear (Fisher) dis- 
criminant JF is derived from the values of sphericity and 
thrust for the event, and the two Legendre moments Lq 
and 1/2 of the energy flow around the i?-candidate thrust 
axis [13 ■ Finally, the polar angle 9b is distributed as 
sin^ 9b for real B decays, while being nearly flat in cos 9b 
for the continuum. 

We select B° —> D+Tr^ and B" —> D^ K+ candidates 
that satisfy Rl > 0.75, and accept 5° -^ D*'^tt^ and 
S° -> D*-K+ candidates with Rl > 0.8. We mea- 
sure the relative efficiency er^ of the Rl selection 
in a copious data sample of decays B^ -^ D^tt^ 
(JJ- -> K+TT-TT-, K°TT-) and B+ -^ D*°n+ (U*" -> 
_D°7, D'^ -^ K^TT^) in which the kinematics is similar to 
that of our signal events, and find that it is consistent 
with Monte Carlo estimates er^ ~ 70%. The fraction of 
continuum background events passing the selection varies 
between 2% and 15%, depending on the mode. 

We identify the signal using the invariant mass m{Ds) 
of Ds candidates and two kinematic variables rriES and 
AE. The first is the beam-energy-substituted mass 
TOES = a/(s/2 -I- pi ■ pB)'^/Ef - p|, where ^/s is the to- 
tal cm. energy, (Ei,pi) is the four-momentum of the 
initial e^e^ system and ps is the B'^ candidate momen- 
tum, both measured in the laboratory frame. The second 
variable is AE = Eg — -\/s/2, where Eg is the B° can- 
didate energy in the cm. frame. For signal events, the 
TO.ES distribution is gaussian centered at the B meson 
mass with a resolution of about 2.5 MeV/c^, and the AE 
distribution has a maximum near zero with a resolution 
of about 17 MeV. The invariant mass m{Ds) has a res- 
olution of (5 — 6) MeV/c^, depending on the D+ decay 
mode. We define a fit region 5.2 < toes < 5.3 GeV/c^, 
|A£;| < 36 MeV, and \m{D s) - m(D s)mG\ < 50 MeV/c^ 
for 5° -^ DfiT" and B'^ — > D^K^ candidates, where 
TO,(Z3s)pDG is the world average Ds mass [lj|. For 
S" -> D*+TT- and S° -^ D*' K+ , we require \m{Ds) - 
to(D,)pdg| < 10 MeV/c2. 

Less than 20% of the selected events in the 
50 ^ D*+Ti- and B° -^ D*- K+ channels (< 4% in 
B^ -^ D^TT^ and 5° -^ D^K~^) contain two or more 
candidates that satisfy the criteria listed above. In such 
events we select a single B^ candidate based on an event 
X^ formed with m{Ds) and Am{D*~^) and their uncer- 
tainties, and the AE variable. Such selection does not 
bias background distributions significantly. 

Four classes of background contribute to the fit region. 
First is the combinatorial background, in which a true 
or fake D)* candidate is combined with a randomly- 
selected pion or kaon. Second, B meson decays such as 
BO^D^*'>+Tr-' , p- with D+- 



^K^'K+ or K*"'K+ can consti- 



tute a background for the B^ -^ Ds' n~ modes if the 
pion in the D decay is misidentified as a kaon {reflec- 
tion background). The reflection background has nearly 
the same thes distribution as the signal but different dis- 
tributions in AE and m[Ds). The corresponding back- 
grounds for the B° -^ DjK+ mode {B°^D-K^*^+) are 
negligible. Third, rare B decays into the same final 
state, such as B°^K'^*'>°K+'!t- or K^*')°K+K- {charm- 
less background), have the same tubs and AE distribu- 
tions as the B° -^ D+tt' or B° -^ D' K+ signal, but 
are nearly flat in m{Ds). The charmless background 
is significant in B'^ — > Dfn^ and B'^ -^ DjK^ decays, 
but is neghgible for B° -^ D*+tt- and B° -^ D*-K+. 
Finally, crossfeed background from misidentification of 
B^ —f Ds' 7r+ events as B^ — > Ds* K^ signal, and vice 
versa, needs to be taken into account. 

We perform a two-dimensional unbinned extended 
maximum-likelihood fit to the toes and m{Ds) distribu- 
tions to extract B{B° -^ D+ti~) and B{B° -^ D~ K+) 
and constrain the contributions from charmless back- 
ground modes. Charmless backgrounds are negligible for 
B" -> D*+TT- and S" -> D*-K+, and we determine the 
BFs of these decays with a one-dimensional fit to the toes 
distribution. For each B decay, we simultaneously fit 
distributions in three Z)+ decay modes, constraining the 
signal BFs to a common value. The likelihood function 
contains the contributions of the signal and the four back- 
ground components discussed above. The combinatorial 
back grou nd is described in to.es by a threshold func- 
tion [1^, dN/dx ex a;\/l-2a;2/sexp [-£, (l - 2x'^ / s)] . 
In m{Ds), the combinatorial background is well described 
by a combination of a first-order polynomial (fake -D^ 
candidates) and a gaussian with (5 — 6) MeV/c^ resolu- 
tion (true D+ candidates). The charmless background is 
parameterized by the signal gaussian shape in tties and 
a first order polynomial in m{Ds). 

For B° -^ D+TT- and B° -^ Dj K+ decays, the fit con- 
strains 14 free parameters: the shape of the combinato- 
rial background ^ (1 parameter for all Df modes), the 
slope of the combinatorial and charmless backgrounds 
in m{Ds) (3 parameters), the fraction of true Df can- 
didates in combinatorial background (3), the number 
of combinatorial background events (3), the number of 
charmless events (3), and the BF of the signal mode 
(1). The signal yields for each Df mode are expressed as 



Nsigi = Ng^Bsig Bi Ei, where N^^ = 230 x 10^ i3, is the 
D+ BF for the mode, Si is the reconstruction efficiency, 
and Bsig is the BF (fit parameter) for the decay. For the 
B° -^ D*+TT~ and B° -^ Dl~K+ decays, 5 free parame- 
ters are determined by the fit: ^ (1 parameter for all Z3+ 
modes) , the number of combinatorial background events 
(3), and the BF of the signal mode (1). The BFs of the 
channels contributing to the reflection background are 
fixed in the fit to the current world average values [lj|, 
and the BFs of the crossfeed backgrounds are determined 
by iterating the fits over each B decay mode. The results 
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FIG. 2: (a,c,e,f) jties projection of the fit with \m{Dt) — 
m(_D+)pDG| < 10 MeV/c^ and (b,d) m{Ds) projection with 
5.275 < ruES < 5.285 GeV for (a,b) B° ^ ZJ+tt", (c,d) 
B" -> D~K+, (e) B° -^ Dl+n-, and (f) B° ^ Dt' K+ . The 
black solid curve corresponds to the full PDF from the com- 
bined fit to all Df decay modes. Individual contributions are 
shown as solid red (signal PDF) , green dashed (combinatorial 
background), and blue dotted (sum of reflection, charmless, 
and crossfeed backgrounds) curves. 



of the fits are shown in Fig.[21and summarized in Tabled 
The systematic errors are dominated by the 13% rela- 



tive uncertainty for B{D^ 



") [l3- The uncertain- 



ties in the relative BFs B{D+^K*°K+)/B{D- -x/jtt ' j 
and B{D+^KlK+)/B{D+^(j)TT+) contribute (5 - 
7)%, depending on the decay channel. Uncertain- 
ties in the selection efficiency are estimated to be 
3% for B° -* D+TT- and B° ^ DjK+, and 7% for 
B° -^ Dl+TT-and B° -^ D*-K+. The uncertainties in 
the refiection and crossfeed backgrounds are below 1% 
for all decay channels. The rest of the systematic errors, 
which include the uncertainties in tracking, photon and 
K'^ reconstruction, charged-kaon identification efficien- 
cies, and variations of the PDF shapes between data and 
Monte Carlo, amount to (6 — 7)%. 

The ratio Phkg = 'Co/^Cmax, where £max is the max- 
imum likelihood value, and Co is the likelihood for a 
fit with the signal contribution set to zero, describes 
the probability of the background to fluctuate to the 
observed number of events. Including systematic un- 
certainties and assuming gaussian-distributed errors, it 
corresponds to the significance of signal observation of 



5 {B° -^ Dfn-), 6 {B° -> L>;+7r- 
and 5 {B° -^ D*' K+) standard 



), 9 (5° - 
deviations. 



D-K^ 
This 



TABLE I: The number of reconstructed candidates (A'raw), the signal yield (A^'sig), computed from the fitted branching fractions, 
combinatorial background (A'comb), and the sum of charmless, reflection, and crossfeed contributions (A'^peak), extracted from 
the likelihood fit. Also given are the reconstruction efficiency (e), the probability (-Pbkg) of the data being consistent with the 
background in the absence of signal, and the measured branching fraction B. The first uncertainty is statistical, and the second 
is systematic. 
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1.32 ±0.27 ±0.15 




D+^KlK+ 


19 


10 ±2 


12 ±4 


1 


9.6 












Dt^(t>TV+ 


197 


32 ±5 


151 ±13 


8±6 


23.4 










B" ^ D-K+ 


Dt->K*'>K+ 


331 


27 ±4 


306 ± 18 


-4 ±6 


17.6 


3 


^0-19 


2.5 ±0.4 ±0.4 
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the first observation of S" -^ Dfir , i?" -^ D*^Tr , and 
B" -^ DI-K+ decays. 

The BF results are collected in Tabled Since the dom- 
inant uncertainty comes from the knowledge of the Df 
BFs, we also report the products B x B{Df -^ (/)7r+). 

The BFs for B^ -^ Ds* K^ are small compared to the 
dominant decays B^ — > D^*'~Tr~^ , implying relatively in- 
significant contributions from the color-suppressed W- 
exchange diagrams. Assuming SU(3) relation, Eq. |^, 
we determine r(i:)7r) = (1.3±0.2(stat)±0.1(syst)) x 10-^ 
and r{D*Tr) = (1.9± 0.2(stat) ± 0.2(syst)) x 10"^^ which 
implies small CP asymmetries in B^^fD'^*'>^TT^ decays. 
These results supersede our previously published mea- 
surements Q- 
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